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Abstract 

Neutrinos can polarize a medium due to their weak interaction. This can manifest itself as an effective 
induced charge of the neutrino. We show how it is related to the Debye screening length in a plasma, 
first using the results of the 1-loop calculation of the neutrino electromagnetic vertex and then in general 
to all orders in e and to leading order in the Fermi coupling. We also discuss how the results are modified 
if the neutrinos have mass ~ either Dirac or Majorana type. 



PACS numbers: 14.60. G 



1 Introduction 



The properties of neutrinos that propagate through a medium have been the subject of great interest in the 
recent Hterature. This has been motivated by the attractive suggestion that the solar neutrino problem 
can be solved by the resonant oscillation mechanism , which hinges on the characteristics of the neutrino 
propagation in a background medium Inspired by its potentially important effects, the neutrino 
interactions in a material environment have also been studied in some detail ^, |6|, |7|, ^, ^, 1^, 11 , ^] . Of 
primary interest along these lines is the study of the electromagnetic interactions of neutrinos in a medium. 
A classic problem in this field is the decay of a plasmon into pairs [^] , which has received considerable 
attention recently ||8|, |ll|, 14]. It has also been pointed out that the rate of decay of a massive neutrino into 
a lighter neutrino and a photon increases tremendously ||9|, |l^ in matter, as a consequence of the fact that 
the GIM mechanism is not operative in a medium with electrons but no muons or taons. 

The study of the electromagnetic interactions of neutrinos in a medium, as well as the results and 
conclusions mentioned above, are based on the 1-loop calculation of the effective electromagnetic vertex of 
the neutrino, which has been performed to the leading order in the Fermi constant ^ using the methods 
of "Quantum Statistical Field Theory''^. The implications of this calculation have been only partially 
explored. In this article, we show that the results of Ref. |p imply that the neutrino acquires a small 
effective charge in a medium. This point was realized by Oraevsky and Semikoz using methods of 
plasma physics before the calculation of the vv^ vertex was performed. On the other hand, our method 
is entirely different, hopefully easier to follow for a particle physicist, and it brings out some interesting 
points which are not easy to see in the method used by Oraevsky and Semikoz Q. In addition, we also 
calculate the induced charge for neutrinos which are massive, distinguishing between the cases of Dirac 
and Major ana masses. 

The plan of the paper is as follows. In Section 2 we establish notation and the definition of the induced 
charge that is used in the rest of the paper, and show why the neutrino can acquire an induced charge in 
the medium but not in the vacuum. Then, using field-theoretic methods, we will show that the neutrino 
electromagnetic vertex is related to the photon self-energy in the medium and, in particular, the neutrino 
induced charge is related to the Debye screening length in a plasma. This relation is derived to 1-loop 
order in Section 3, and then in Section 4 we show that it is valid to all orders in e (and to first order in 
the Fermi coupling). In Section 5 we use this relation to find an expression for the induced charge of the 
neutrino, and in Section 6 we show that the magnitude of the induced charge depends on whether the 
neutrinos are massless or massive and, in the latter case, on whether the mass is Dirac or Majorana type. 
Finally, using the results of Ref. [0], we will estimate the induced charge in some particular backgrounds. 



2 Electromagnetic vertex of the neutrino and the definition of the in- 
duced charge 

We begin by understanding the reason why the neutrino can acquire an induced charge in the medium 
but not in the vacuum. The off-shell electromagnetic vertex function F;^ is defined in such a way that, for 
on-shell neutrinos, the vv^ amplitude is given by 

M = -iu{k')Txu{k)A^{q) , (2.1) 

where 



q = k-k' (2.2) 

^More often, this is called the "Finite Temperature Field Theory", but the name is misleading because the methods are 
also applicable at zero temperature with a finite density of particles in the background. 
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is the momentum carried by the photon. In general, Tx depends on k and k' or, equivalently, on k and q. For 
neutrinos in a medium T\ depends also on the parameters characterizing the medium. For homogeneous 
and isotropic media, to which we will restrict ourselves, there is only one such parameter, viz. the velocity 
4- vector of the background medium v^. 

There are two important consequences of the fact that the external (neutrino) lines in the Feynman 
diagram for the viyj amplitude are neutral. Firstly, Tx satisfies 

qxT'' = (2.3) 



for all values of q. It is important to realize that for neutrinos Eq. (|2.3| ) holds for arbitrary values of q, and 
not just when k and k' are on shell. If the fermion lines in the diagram were to correspond to a charged 
particle (e.g., the electron), then the analogous relation in that case involves terms in the right-hand side 
involving the inverse propagators corresponding to the external fermions. In that case Eq. ( p.3| ) does not 
hold for arbitratry values of q, but only when both k and k' are on shell. 

The other important consequence of the fact that the external lines are neutral is that Tx is well defined 
in the limit q^ — > 0. The reason is that the photon vertex must be connected to a pair of internal lines of 
the diagram. If one of these lines is assigned the loop momentum p which is integrated over, the other line 
will carry a momentum pig. The propagator of this second line will involve the factor 

(2.4) 



{p lb q)"^ — m? q"^ ziz2p ■ q + {p"^ — m?) ' 

where m is the mass of the internal line. However, since p^ ^ w? for any internal line, no singularity is 
produced for q^^ 0. 

From the two properties of F;^ just discussed, we obtain 

FA(g'^ = 0) = 0, (2.5) 

which implies that the particle associated with the external line does not acquire a charge in any order 
of perturbation theory. To see this explicitly, we will consider the matrix element of the charge operator 
between two neutrino states with momenta: 

k^ = {E,k), k'^ = {E,k'). (2.6) 

We use states with the same energy, because then 

q = k-k' = {0,q) (2.7) 

with q = k — k', which corresponds to the static limit. Denoting by p{x) the charge density operator, the 
effective charge is defined by the equation 

ees{k'\k) = I d^x{k'\p{x)\k) 

= i27rf 6(-^\q){k'\pmk) (2.8) 

On the other hand, 

{k'\p{0)\k)=u{k')To{0,q)u{k), (2.9) 

where the notation Tx{q^,q) has been used to indicate explicitly that we are considering the dependence 
of F separately on the frequency and wavelength of the photon. Thus we obtain 

eeff = ^u{k)To{0, q^ 0)u{k) , (2.10) 
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which is the basic equation to interpret our results, but it can be cast in an elegant way. Introducing the 
spinor projection matrix 

Sik)=u{k)<S)u{k) = ^{l + Xj5W, (2.11) 

where the last step is valid for massless Weyl spinors with A = ±1 being the helicity, we can rewrite Eq. 
( |2l0| ) as 

eeff = ^tr[ro(0,g^O)5(fc)], (2.12) 

= ^ tr [To{0,q^ 0)(1 + A75)^] , (2.13) 

where again the second step is valid for massless Weyl spinors. 

Since Tx has a well defined limit as g'' — > 0, we can make a Taylor expansion around = 0: 

To = Go + g°Gi + 0((?o') 

f = #o + g°^i + 0('?°'), (2.14) 



where all the coefficients are independent of q^. Then from Eq. ( |2.3| ) it is easy to deduce that, in the limit 

To = q-Hi + 0{q^) 
f = q'^Hi+0{q'''). (2.15) 

Since Hi has a well defined limit as g — > 0, it follows that Eq = in this limit, and from Eq. ( p.l3| ) we see 
that e^"^^ is zero in the vacuum. 

In the medium, Eq. ( |2.3D continues to hold for any value of q, and the relations in Eq. ( 2.15| ) are also 



valid, but Hi no longer has a well defined limit as g ^ 0. There are several ways to understand why this 
is so. One of them is to notice that in the case of the medium, some of the internal lines to which the 
photon is attached are on shell because they correspond to particles that are in the background. Thus, the 
singularities that are avoided in the case of the vacuum because the photon is attached only to internal 
off-shell lines, reappear here. Therefore, nothing prevents Hi to develop a singularity as g ^ of the form 

Hi = (constant) • rSr (2-16) 
q^ 

In such a case, the constant appearing in this equation is the value of Eq in the limit g'^ = 0, g ^ 0. Thus, 
from the definition of the effective charge in Eq. ( ^.13 ), it follows that e^g is non- vanishing. 



In what follows, we will show that this is exactly what happens, first explicitly by using the 1-loop 
calculation of the neutrino electromagnetic vertex and then by a general field-theoretic argument, which 
extends the 1-loop result to all orders in e. This will be done by showing that Tx is related to the photon 
self-energy TTxp{q)- Eurther, in the limit that we are considering, T^xpiq^ = 0, g — > 0) is related to the Debye 
screening length, which allows us to establish the relation between the latter quantity and the neutrino 
induced charge. 



3 One-loop result 

We are interested in the regime where the neutrino momenta are small compared to the masses of the 
W and Z bosons. Therefore, we can neglect the momentum dependence in the W and Z propagators. 
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which is justified if we are performing a calculation to the leading order in the Fermi constant Gp- In this 
approximation, the diagrams contributing to the electromagnetic vertex then appear at the 1-loop level, 
and are shown in Fig. ||. Since the momentum dependence of the weak gauge bosons are neglected, these 
two diagrams can be represented by the single diagram of Fig. |2| with a four-fermion vertex. Let us denote 
the 4-fermion interaction by 

^(weak) ^ ^j^^^^ ^ ^^^^^j ^3^^^ 

where L = ^(1 — 75) is the projection operator for left chirality, and / stands for the electron field. We 
can then write the amplitude of Fig. Q as 

-^Fa = {-ie){-iGF^2){-l)YL X 

/ iSF{p)^p{A + ^75) iSF{p -q)], (3.2) 

where iSpip) denotes the propagator of the background particles with momentum p, and e is the charge 
of the electron. In complicated systems, this propagator may be complicated, and the integration over the 
momentum p may have unusual measure as well, but we do not need these explicitly for what follows. 

Now the interesting observation is that the contribution from A in Eq. ( |3.2D is intimately related to 
the vacuum polarization of the photon which, at 1-loop, arises from the diagram in Fig. and is given by 



i-^Xpiq) = (-^e)^(-l) / 77^^ tr ['yxiSFip)'ypiSFip - q)] ■ (3.3) 



d'^p 

W) 

Therefore, Eq. (|3.2[) can be written in the form 



.^Y{l-l'){AiiXp + Birl,). (3.4) 



where we have defined 



^ttap = (-«e)^(-l) j -^i^^"^ h\iSF{php75iSF{p - q)] ■ (3.5) 

The term proportional to tt^^ does not contribute to Fq. The reason is that the trace contains a factor 
of 75, and therefore can be non-zero only if there are at least four other 7-matrices present. Since the 
terms in a fermion propagator have at most one 7-matrix, the trace involves tr (757A7P 70,7/3) = 'iie\paf3- 
After the p integration this can yield only a term proportional to expapq'^v^ , which does not contribute to 
the zeroth component of F^. 

4 Generalization to higher orders 

The above result, based on the 1-loop calculation of the photon self-energy and the neutrino vertex function, 
can be generalized as follows. The photon self-energy is defined in general by 



= {-ief j d''xe-'^-^iTU^-^^3xmp{x)), (4.1) 
where jx is the electron current density 

jx=llxf, (4.2) 
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and 

U^'^^ =exp(^il d^'zClX^^ (4.3) 

with 

ctf = -enA\ (4.4) 



On the other hand, the neutrino vertex function of Eq. ( p.l| ) is defined by 

r^{k,k') = ej d''xdSe-'^--e^^'y(TexY>{i J 4^^^) z.(2/)jA(0)F(x)^ , (4.5) 

where 

^ (total) _ ^(weak) ^(em) / , „x 

''-int — ''-int ''-int • \^-^) 

The subscript a in Eq. ( |4.5| ) is used to indicate that Tx is obtained from the above formula by amputating 
the propagators corresponding to the external neutrino lines. It is convenient to rewrite c[^^^^\ given in 
Eq. ( |3.1] ), in the form 

where 

jl^Jlxl'f- (4.8) 

To first order in Gp, 

rx = ej d^d^yd^ze-^'^-^e^'^'-y {TU^'^"^^v{y)jx{mx)iCtr''\z)) ^ (4.9) 
which, using Eq. ( [d.T] ) and amputating the external neutrino lines, reduces to 

Tx = -ieGFV2AYL X 

j d^ze-'^i-^ (Ti7(<=-)iA(0)jp(z)) + j d^ze-'^i-' (t C/(^-) jA(0)i,^(z))) . (4.10) 



Defining 



i-ief j d^xe-'^-^ {TU^^^^h{^)f;{x)) , (4.11) 



Ta = -:^7''(1 - + Bnl) . (4.12) 



we finally obtain the relation 



This expression is the same as Eq. ( |3.4D , except that now it is clear that it is valid in all orders of 
the electromagnetic interactions. Since vr^^ is a pseudotensor that depends only on q and v, it must be 
proportional to expapq^v^, and therefore it does not contribute to Fq in the rest frame of the medium. If 
one includes the effects of strong interaction as well, the various occurences of C^^^ should be replaced by 
£^^^ + C'^'°''^\ but this would not change Eq. (p^) . 
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5 Expression for the induced charge 

The relation between the induced charge and the Debye screening length is obtained as follows. As already 
argued, the zeroth-component of Tx is given by 

ro(0,g - 0) = -^7^[l - 7']vrop(0,g - 0) . (5.1) 
The most general form of nxp is [ p!5| ] 

vTAp = t^tRxp + t^lQxp + T^pPxp , (5.2) 

where 

R\p = 9Xp — Q\p , (5.3) 

Qap = (5.4) 



P\p = -^^Xpapq'^v'^ , (5.5) 



with 



9\p = gxp ^ (5.6) 

vx = gxpv'' (5.7) 

Q ^ ^(g.t;)2-g2. (5.8) 

In the rest frame of the medium, with = (1,0,0,0), the above definitions give 

Qxp{0,q) = VxVp, (5.9) 

Roo{0,q)=Rio{0,q)=Ro,{0,q) = 0, (5.10) 

PooiO,q)=PioiO,q)=PoiiO,q) = 0, (5.11) 

where we have indicated explicitly the fact that we are evaluating R, Q and P in the static limit, = 0. 
From Eq. (5.2) we then obtain 

7roo(0,g) = 7rL{0,q), (5.12) 

TToMq) = ^^^o{0,q)=0. (5.13) 

Eq. (|5.1]) then yields 

ro(0,g ^ 0) = - (^^) 7o(l - 75)vrL(0,g - 0) . (5.14) 
Using the definition of the Debye screening length, 

r^^ = TTLiO,q^O), (5.15) 



we finally obtain for the induced neutrino charge by the use of Eq. (2.13): 

where, as stated before, A is the helicity of the neutrino. Thus, it is clear that only the left-handed neutrinos 
have an induced charge. The induced charge for the right handed neutrinos vanish since they have no weak 
interactions. If they interact via some other weaker interaction, then of course they also acquire an induced 
charge, but the magnitude of that will be further suppressed. Also, note that Cj^^ oc e, since r|) oc e ^ 
which follows from Eq. (|5.15|). 



7 



6 Generalization to massive neutrinos 



The generahzation of our previous results to massive neutrinos is straightforward, ahhough there are 
several important differences. The effective charge is defined by Eq. ( p. 12 ), but the expression for the 



spinor projection operator S{k), as well as the expressison for Fq, differ from the previous case. The spinor 
projection operator S{k) is given, for massive particles, by 

S{k) = u{k) ® u{k) = ^(^ + m)(l + A7V) , (6.1) 

where is the spin polarization vector which, for helicity states, is given by 

s^^ = -(\k\,Ek/\k\) . (6.2) 

Therefore, although Eq. ( p.l2| ) remains valid for massive neutrinos, Eq. ( 2.13| ) does not. To proceed, we 
consider the cases of Dirac and Major ana neutrinos separately. 

6.1 Dirac case 



For Dirac neutrinos, F;^ remains to be given by Eq. (4.12) and the effective charge by Eq. ( [2.121 ). The 



formula for the effective charge of Dirac neutrinos is now obtained by substituting Eq. (6.1) in Eq. ( 2.12| ), 
yielding 



(.15) ^ _GfA_ ( _ 

V2erl [ E , 



(6.3) 



Specializing this formula to the case of massless neutrinos we recover Eq. ( ^.16 ), as it should be. In the 



opposite limit of non-relativistic neutrinos, we see that both helicity states have the same effective charge, 
which is equal to half the value of the charge of the left-handed neutrino in the massless case. 

6.2 Major ana case 

For (massive) Majorana neutrinos we again have to use the spinor projection operator appropriate for 
massive particles, given in Eq. ([6.1[ ). However, the formula for Fa is modified as follows. Going back to 
Eq. (4.9), it is important to recognize that each one of the neutrino field operators y{y) and V{y) can be 



contracted with either one of the same field operators that come from the factor C^^^'^^\ The reason is 
that for a Majorana neutrino the field operator is self-conjugate and therefore ^{y) can be contracted with 
not only v{y) but also with itself. The net result of adding these two possible contractions is that the 
expression for Tx given in Eq. ( ^1^ is replaced by |l§ 



= -:^(-27V)(-4vrv + 37.%) . (6.4) 



Substituting this formula and the projection operator given above for the massive case, intoEq. ( ^l2[ ) we 
then obtain the effective charge for Majorana neutrinos 



^(.M) _ V2GfA \\k\ 
'''''^ ~ erl E~ 



(6.5) 



We notice the following features: (i) the positive and negative helicity states have opposite effective charge; 
(a) in the masless limit the result for the negative helicity state is the same one obtained in the massless 
Dirac case and in the Weyl case, while the positive helicity state has the opposite value of the charge. This 
is not surprising since the right-handed component of the Majorana neutrino is just the CPT conjugate of 
the left-handed one. 
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7 Numerical estimates 



In order to obtain numerical estimates, we note that since \k\/E < 1, the magnitude of the induced charge 
is maximum if the neutrino is massless. Thus, in this section, we use Eq. ( 5.16| ), which is vahd for massless 



neutrinos. We see that we need to know two things in order to obtain a numerical estimate for the induced 
charge of the neutrino, viz., A and r^. The first is easy, and is given by the standard model of electroweak 
interactions. In fact, one obtains 

1 2sin^^;^ + i for., 

1 2 sm &w — ^ tor f^, Ur . 

For massless i^e, our result exactly reproduces the results of Oraevsky and Semikoz [^]. Our formulas can 
be used to obtain the induced charges for the and ly-r as well, even if they are massive particles. 

One curious thing to note is the fact that the induced charge of the Ve is different from that of or 
Uj- since the i^e's interact with the electrons of the medium via both charge and neutral currents. Thus, if 
neutrinos mix, when they oscillate during their propagation through a medium, the induced charges also 
oscillate. Using sin^ 9\y = 0.23, we see that 



-24. (7.2) 



ind 



It therefore seems that this oscillation of charges should be a fantastic phenomenon, judging by the fact 
that the ratio of the induced charges is large, and also is negative. However, that is not the case because 
the magnitude of the induced charge appears to be extremely small. In fact, putting numbers in Eq. ( 5.16| ), 
one obtains 

4"5' = -2.10-xfiS^)'. p.3) 



To proceed, we need the value of ro for the background. This can be obtained either from the results of 
Ref. 1^] or from standard texts on plasma physics. For a background consisting of non-relativistic electrons 
at temperature T, the Debye radius is given by 

rl = ^, (7.4) 

where is the electron number density^. In order for the induced neutrino charge to be detectable in 
experiments, the values of T and rie must be such that the resulting Debye radius is small enough. While 
this is not the case for any known plasma, the methods presented here may prove useful in applications to 
similar situations where more exciting results may be obtained. 

Note added in proof : After this work was submitted for publication, we were made aware of a paper 
by Altherr and Kainulainen [0] where the one-loop electromagnetic vertex of neutrinos has been calculated 
in a medium. The calculation agrees with that of Ref. Q. These authors specifically noted that an induced 
charge appears in the medium. No effort was made to make contact with the Debye radius. 

Acknowledgements : The work of PBP was supported by a grant from the Department of Energy. 



^Formulas appearing in books on Plasma Physics usually have a factor 47r in the denominator on the right hand side since 
their definition of electric charge is different. 
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Figure 1: 1-loop diagrams for the effective electromagnetic vertex of neutrinos which contribute in the limit 
that the 4-momenta of the W and the Z lines are neglected. 




Figure 3: 1-loop diagram for the vacuum polarization of the photon. 
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